The search for 0-N bound state from subthresold production of (p meson 
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The subthreshold photoproduction of <j) meson from heavy nuclear targets has been suggested 
as a candidate to search for the i^-N bound state, a quantum chromodynamics (QCD) molecular 
state. In this paper, we present detailed Monte Carlo studies to demonstrate the feasibility of this 
technique. Further, we show that proton induced subthreshold production of (j> meson from heavy 
nuclear targets is also suitable for such a search. 
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S. J. Brodsky, I. A. Schmidt and G. F. de Teramond 
suggested that the QCD van der Waals interaction, me- 
diated by multi-gluon exchanges, is dominant when the 
two interacting color singlet hadrons have no common 
quarks. Luke, Manohar and Savage 0] predicted that the 
QCD van der Waals interaction is enhanced at low veloc- 
ity. This finding supports the prediction that a nuclear- 
bound quarkonium can be produced in charm produc- 
tion reactions at threshold. Brodsky, Schmidt, and de 
Teramond [Ij] investigated the nuclear-bound quarko- 
nium state using a non-relativistic Yukawa type attrac- 
tive potential V^gg^^ = —ae~^'^/r characterizing the 
QCD van der Waals interaction. They determined the 
a and /i constants using the phenomenological model 
of high-energy Pomeron interactions Using a vari- 
ational wave function ^'(r) = (7'^/7r)^/^e~''''', they pre- 
dicted bound states of rjc with '^He and heavier nuclei. 
This prediction was later confirmed by Wasson [4] using 
a more realistic V^gg^^ potential taking into account the 
nucleon distribution inside the nucleus. 

Similarly, one expects the attractive QCD van der 
Waals force dominates the 0-N interaction since the 
(j) meson is almost a pure ss state. Using the varia- 
tional method and following Ref. [jj to assume Vs^ssYn = 
—ae^^/r, Gao, Lee and Marinov 'gj find that a bound 
state of (/)-N is possible with a = 1.25 and /i = 0.6 GeV. 
The binding energy obtained is 1.8 MeV. Their results 
should be compared with a = 0.6 and /i = 0.6 GeV de- 
termined in [l|] for the cc quarkonium. The interaction is 
expected to be enhanced by [mc/msf' , i.e., qq separation 
cubed, from cc to ss. Since the radius of the <j) meson is 
0.4 fm twice the radius of the J/^! meson, a — 1.25 is 
a rather conservative coupling constant to use for the (jy 
N interaction. Also, the interaction is expected to have 
longer range for the 0-N system than that of the cc-N 
interaction. Thus, /i = 0.6 GeV used in the variational 
approach described above is also conservative for the cjy-'H 
interaction. 

Recently, the 0-N bound state has been studied by 
Huang, Zhang and Yu 7] using a chiral SU(3) quark 
model and the extended chiral SU(3) quark model solv- 
ing the Resonating Group Method (RGM) equation. The 



model parameters used by the authors in this work pro- 
vided good description of baryon bound states, deuteron 
binding energy and NN scattering phases shifts in their 
previous work d, Q . A 0-N quasi-bound state with sev- 
eral MeV of binding energy was predicted by the afore- 
mentioned extended chiral quark model plus channel cou- 
pling effect. 

Such a 0-N bound state could be formed Q inside a 
nucleus at the quasifree subthreshold photoproduction 
kinematics where the attractive force is expected to be 
enhanced. Recently, "subthreshold" (f) meson photopro- 
duction has been observed 10] for the first time from 
a deuterium target at Jefferson Lab. The experimental 
search for such a bound state would be a triple coinci- 
dence detection of kinematically correlated , , and 
proton in the final state from subthreshold production of 
(j) meson from heavy nuclear targets. To identify clearly 
the signal of a (/)-N bound state, one needs to understand 
the background contribution to the pK~^K~ final state 
carefully. The dominant background contribution to the 
signal of interest are the quasifree subthreshold produc- 
tion of (j) meson from a bound proton and the subsequent 
decay of (j> into K'^ K" without the formation of a (/)-N 
bound state, and the direct production of from 
a bound proton inside the nucleus. Recently, we car- 
ried out a detailed Monte Carlo simulation of these pro- 
cesses. The Monte Carlo study shows that cuts on the 
momentum correlation between proton and , and on 
the invariant mass of the , and proton system can 
clearly separate the signal of the decay of a 0-N bound 
state from the backgrounds. Therefore, one can identify 
a bound 0-N state experimentally using the aforemen- 
tioned triple coincidence experimental technique. In this 
paper, we present our results from the Monte Carlo stud- 
ies. 

The kinematics for all three channels considered in the 
Monte Carlo study for the detection of a (/)-N bound 
state, using ^^C, ^^Fe, ^'^Cu and ^^^Au target nuclei 
and the CLAS detector [llj at Jefferson Lab, follow a 
set of common parameters and assumptions. First, the 
energy of the photons in the region of interest is dis- 
tributed uniformly from 1.65 to 1.75 GeV. This energy 



2 



range is roughly 80-180 MeV higher than that of a sim- 
ulation in the case of an ideal detector for the triple co- 
incidence detection of proton, K^,K^. The target nu- 
clei are assumed to be initially at rest in the lab frame. 
For each event, the Fermi momentum and the missing 
energy for the bound proton are weighted by the nu- 
clear spectral function. The nuclear spectral functions 
for the ^^Fe, ^'^Cu, and ^^^Au nuclei used in our simu- 
lations are based on the Green's function Monte Carlo 
calculation and from de-radiating data for the ^^C 
target A Breit-Wigner distribution with a mean 

value of 1.019413 GeV/c^ and a width 0.00443 GeV/c^ 
is used to model the mass of the </> meson that is pro- 
duced. Furthermore, the mass for each event is bound 
between 0.819413 to 1.219413 GeV/c^ to avoid physically 
unreasonable masses. 

Furthermore, the CLAS detector is modeled by apply- 
ing cut-offs to the zenith lab angles that can be detected 
for the final state of each channel. From an analysis of the 
CLAS glO data for the final state oipK^K^ from me- 
son photoproduction [l3| , the minimum zenith lab angles 
for and K~ are 5 and 20 degrees, respectively. Ad- 
ditionally, the resolution of the detector is incorporated 
into the simulation by weighing events using a Gaussian 
type distribution following the procedure developed for 
the Jefferson Lab Hall C proposal on electroproduction 
of (f) using Hall C SOS and HKS magnets employing the 
missing mass technique [l5| . Realistic numbers for the 
CLAS detector resolutions are used in our simulation. 

The first channel considered is the <f> meson produc- 
tion from a bound proton with the subsequent decay of 
(f) into K^K^ without the formation of a (/)-N bound 
state. The events simulated are weighted by the mass 
of the (j) meson, and the Fermi momentum and missing 
energy of the bound proton. Before computing the kine- 
matics of 7 -t- "p" p + (j), the energies of 7 and "p" 
are checked to ensure they are sufficient to produce a 
(j) meson; events with insufficient energy are discarded. 
Given that no bound state is formed, the kinematics of 
the (j) meson decay into K'^K~ are calculated. At this 
point in the simulation, events for which the zenith lab 
angle of or K~ below the CLAS detector cut-offs 
are removed. Before simulating the detector's resolution, 
the calculations performed are tested in two ways: recon- 
structing the mass of the (j) meson from the energy and 
linear momentum of the pK~^K~ final state; and recon- 
structing Fermi momentum of the bound proton from the 
4-vector of the initial and final state particles. Finally, 
the detector resolution is simulated for the pK^K~ final 
state. 

In the case of the second background channel, the di- 
rect production of K~ from a bound proton, the sim- 
ulation's structure is essentially the same as the struc- 
ture of the first channel considered. An important dis- 
tinction for the second background channel is the di- 
rect production of the three-particle final state from the 



7 -I- "p" ^ p -\- K'^ -\- K~ process, which is simulated by 
the following sequence of steps. For computational pur- 
poses, it is assumed that the direct production of K^K^ 
can be simulated by 7 "p" {AB) + C ^ A + B + C, 
where B and C are combinations of proton, and 
. The intermediate step of the simulation, {AB) + C, 
has no physical significance, and merely serves as a tool 
for kinematic calculations. Given that are assumed 
to be kinematically indistinguishable in the simulation, 
there are only three kinematically distinct combinations, 
each of which is assumed to have an equal probability of 
occurring. Thus, each event is weighted by the uniform 
distribution of the different sequence of the same three- 
particle final state. Similarly to the first channel, the 
Fermi momentum of the bound proton is reconstructed 
to ensure all calculations are carried out correctly. 

The last channel being studied is exactly the same as 
the first, with the exception of the formation of a (/)-N 
bound state in the nuclear medium. For the purposes 
of the simulation and the detection setup, the nucleon 
to which the 4> meson binds is modeled as a proton. It 
is assumed that QCD van der Waals forces lead to the 
formation of a 0-N between the photo-produced cj) meson 
and a bound proton if the cj) particle has a linear mo- 
mentum less than 400 MeV/c. The momentum of the 
bound proton is assumed to be the same as that of the 
(j) when considering the formation of the 0-N state. The 
0-N exotic state is assumed to have a binding energy 2.5 
MeV [H]. When considering the decay of 0-N into the 
pK^K^ final state, the kinematic methodology used for 
the direct production of K^K~ from bound proton is 
reused. Once again, the simulation is checked by recon- 
structing the mass of the cj) meson and the Fermi mo- 
mentum of the bound proton. Fig. 1 shows the invariant 
mass distribution of the p,K^,K^ system for all these 
three channels, and the histogram for each of the chan- 
nels has been normalized to unity. In the case of the (/)-N 
bound state, the prominent peak corresponds to the pro- 
ton from the decay of the (/)-N bound state is detected 
with the second bump corresponding to the recoil pro- 
ton(s) from the subthreshold <f) production process being 
detected. 

The search of the presence of a 0-N bound state is 
based on the triple coincidence detection of the kinemat- 
ically correlated , K~, and the proton in the final 
state. Such a signature is clearly identified when com- 
paring the absolute value of the linear momentum of the 
and the scattered protons. We are able to segregate 
the events corresponding to the channel with the (f)-N 
bound state by applying one graphical cut, and another 
cut on the invariant mass of the pK^K^ final state. The 
graphical cut involves removing all events for which the 
square sum of the linear momentum of or K'^ and 
the proton is greater than (0.300 GeV/c)^. The second 
cut consists of only considering events where the invari- 
ant mass of the pK~^K~ final state falls in the region 
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FIG. 1: Monte Carlo simulation of the invariant mass distribution of the proton, K~ and K system for the following three 
processes from a copper target: 0-N bound state (red); quasi- free subthreshold production without the formation of a bound 
state (green); direct production oi the pK'^ K~ process (blue). The left and the right panel show the photo- and proton induced 
subthreshold production, respectively. 



1.0 

_0.: 
u 

> 

CD 

a 0.6 - 



I Meson Photoproduction 

Without Bound State 

• Direct Production 

• With Bound State 



0.3 Q.4 0.5 0.6 
Pk-I (GeV/c) 



1 Meson Proton Induced Production 



'01.5- 
> 

CD 

o - 
t~5;i.o- 




3.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 
Pk-I (GeV/c) 



FIG. 2: Monte Carlo simulation of the proton momentum versus K~ momentum for the three processes from a copper target 
described in the text. To improve the visibility of the figure, we plot 25,000 successful events for each channel only, while our 
overall results are based on analyzing 45,000 successful events for each channel. 



of 1.955 ± 0.010 GeV/c^. The results from the simula- 
tion using ^^C, ^^Fe, ^^Cu and ^^^Au target nuclei have 
demonstrated that the channel for 0-N is completely seg- 
regated from the two other background channels. Only 
about 11% of events for the channel containing the (p-N 
bound state are eliminated due to these two cuts with 
no additional cuts from the detectors. An example is 
shown for the case of the Copper nucleus before (Fig. 2) 
and after (Fig. 3) the momentum and invariant mass cuts 
discussed in the text. 

To increase the feasibility of experimental detection of 
a 0-N bound state a Monte Carlos analysis similar to that 
of subthreshold cf) meson photoproduction is carried out 
for subthreshold proton induced (p production. The over- 
all conclusions from the proton induced production are 
the same as those for photoproduction, namely that the 
events containing 0-N bound states can be completely 
separated from background channels. The three chan- 
nels considered for the proton induced (f> production are 
equivalent to those considered in photoproduction, with 



the distinction of the different incident beam. Once again 
the simulation is conducted in the subthreshold region, 
to minimize the momentum mismatch between the pro- 
duced 4> meson and the bound nucleons, and to obtain a 
clear kinematic signature for the presences of i^-N bound 
state when comparing the final state momentum value of 
the final products. The incident proton's total energy for 
the simulated events is weighed by a uniform distribution 
in the subthreshold range of 3.30 - 3.50 GeV, where the 
threshold energy is 3.53 GeV. 

Even though the analysis for proton induced produc- 
tion is similar to that of photoproduction; there are some 
important kinematic distinctions that arise from having 
to consider an additional proton during the various stages 
of the simulation. When simulating the channels where 
a (j) meson is produced from a bound proton the method- 
ology used for the direct production of K^K~ in photo- 
production is reused noting that the different kinematic 
combinations now involve two protons and one meson. 
For the channel corresponding to the direct production of 
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FIG. 3: Same as Fig. 2 except for the momentum and invariant mass cuts applied (see text). 



K^K we need to consider the four-particle final state: 
p + "p" = p + p + + K~ . This process is simulated 
by p + "p" ^ [AB] + {CD) ^ A + B + C + D, where 
A,B,C, and D are combinations of p^p,K^^ and . 
Following the assumption that are kinematically in- 
distinguishable we observed only two dynamically differ- 
ent sequences, which are attributed an equal probability 
of occurring. Each event is weighted by uniform distri- 
bution of the two different sequences. Another impor- 
tant difference of the proton induced <f) production when 
compared to photoproduction is that for the channel as- 
suming the presence of a i^-N bound state the condition 
for such a bound state occurring is that the linear mo- 
mentum of the (f) meson is less 500 MeV/c (as oppose 
to the 400 MeV/c used for photoproduction). This in- 
crease of the (j) meson upper momentum limit is necessary 
due to the overall higher momentum distribution of the 
(j) meson produced. Although high momentum protons 
inside nuclei are suppressed, one still expects more (/)-N 
bound state events from proton induced reaction due to 
the nature of the strong interaction. Such an experiment 
can be carried out in the future at the CSR facility in 
Lanzhou, China where one can design a new detection 
system particularly suitable for this search. 

The presence of a (/)-N bound state is observed when 
comparing the final state momentum of the protons and 
kaons using two cuts to isolate the channel with the (/)-N 
bound state in the proton induced production. The first 
cut removes all the events for which the squared sum 
of the final state momentum of K~ or and proton is 
greater than (0.450 GeV/c)^. The second cut is the same 
as that for the photoproduction case which is on the in- 
variant mass of the final state pK~K'^ . These two cuts 
yield a clear and uncontaminated kinematic signature for 
the presence of a <?!>-N bound state. The percent on events 
from the channel containing the (^-N bound state elimi- 
nated due to these two kinematic cuts is about 11% with 
no additional cuts from the detectors. 

In the two cases presented in this work, the kinematic 
cuts described are able to separate the signal cleanly 
from the main physics backgrounds, therefore it is not 



necessary to carry out the simulations with cross-section 
weighting. To carry out such a search at Jefferson Lab, 
the feasibihty of using the CLAS BONUS ^ type recoil 
detector for detecting low-momentum charged particles 
will be studied in the near future. 

In summary, we carried out detailed Monte Carlo stud- 
ies of subthreshold (j) meson production from photo- and 
proton induced reactions to demonstrate the feasibility 
of experimental search for 0-N bound state from heavy 
nuclear targets at Jefferson Lab in U.S. and the CSR 
facility in Lanzhou, China. This work is supported in 
part by the U.S. Department of Energy under contract 
number DE-FG02-03ER41231. H.G. also acknowledges 
the support of the National Science Foundation of China 
through an overseas young scholar collaborative award. 
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